Dependence of the evolution of the cavity radiation of a coherently pumped correlated 
emission laser on dephasing and phase fluctuation 
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Analysis of the dynamics of the cavity radiation of a coherently pumped correlated emission laser 
is presented. The phase fluctuation and dephasing are found to affect the time evolution of the 
two-mode squeezing and intensity of the cavity radiation significantly. The intensity and degree of 
the two-mode squeezing increase at early stages of the process with time, but this trend changes 
rapidly afterwards. It is also shown that they increase with phase fluctuation and dephasing in the 
strong driving limit, however the situation appears to be opposite in the weak driving limit. This 
essentially suggests that the phase fluctuation and dephasing weaken the coherence induced by a 
strong driving mechanism so that the spontaneous emission gets a chance. The other important 
aspect of the phase fluctuation, in this regard, is the relaxation of the time at which the maximum 
squeezing is manifested as well as the time in which the radiation remains in a squeezed state. 

PACS numbers: 42.50.-p, 42.50.Ar, 42.50.Gy, 42.50.Lc 



I. INTRODUCTION 

In a two-photon three-level laser, the nonclassical fea- 
tures of the radiation are attributed to the atomic coher- 
ence induced by preparing the atoms initially in a coher- 
ent superposition of the energy levels between which a 
direct spontaneous transition is electric dipole forbidden 
[1|4J| or by pumping the same with external radiation [5j- 
KLG ] . The roles of these mechanisms on the generation of 
light with nonclassical features have been extensively in- 
vestigated elsewhere fll-[TTj]. The earlier studies unequiv- 
ocally assert that the various schemes of this laser can 
be a potential candidate for the realization of a strong 
bright entangled light. However, a number of naturally 
occurring phenomena that significantly alter the degree 
of nonclassicality were ignored so far for various reasons. 

To begin with, as recently discussed, the present prac- 
tical capability cannot provide a foolproof for a perfect 
preparation of the atomic coherent superposition [l2| - 
fl5j . In other words, the actually induceable coherent 
superposition is less than what is assumed due to loses 
arising from faulty setups. Even though there is a pro- 
posal acclaiming that the phase fluctuations that emerge 
in connection with the faulty preparation might enhance 
nonclassicality via creating indirect spontaneous emission 
roots [IH , a detailed analysis is still missing to attest the 
claim. The other noteworthy mechanism that has a po- 
tential for destroying the coherence is the dephasing that 
originates from loses of atomic coherence through vari- 
ous quantum processes. In light of this, in a case when 
it is possible to lock the phase, it has been argued that 
there is a chance for enhancing the nonclassicality [l6j . 
When there is no external intervention, it is more likely 
that the atomic coherence decays faster than the atomic 
decay rates. 

Moreover, the study of the time development of the 
two-mode squeezing and entanglement has received a 
considerable attention in recent years 



has been generally observed that the two-mode squeezing 
and entanglement increase at early stages of the process 
and later found to decrease with time [H, [l3|-[l5j . Partic- 
ularly, in the pumped correlated emission laser when the 
atoms are initially prepared to be perfectly in the lower 
energy level and the external pumping radiation induces 
phase fluctuations, the emerging phase fluctuations turn 
out to decrease the entanglement for all cases [TtJ • Be- 
sides, the effect of the phase fluctuations is found to be 
prominent when the cavity is assumed to contain photons 
initially [l3| . Investigation of a similar nature with no 
external driving radiation when the atoms are partially 
prepared indicates that the phase fluctuations induced 
due to imperfect preparation decrease entanglement and 
the time for which the system remains entangled [Til ]. 

With this background, it appears natural to raise some 
basic questions. Is there a chance for the phase fluctua- 
tions to enhance the entanglement? Assuming that it is 
technically possible, what happens if the atoms are exter- 
nally driven by a strong coherent light? Can the degree 
of entanglement increase if the coherent superposition 
decays with larger rate than the spontaneous atomic de- 
cay rate? The fundamental issue behind these questions 
is the suspicion emanating from the competing nature 
of the various phenomena in relation to the coherence 
responsible for the manifestation of the nonclassical cor- 
relations. To answer these and other similar questions, a 
comprehensive approach in which these competing effects 
are brought together is required. 

Hence, in this contribution, coherently pumped corre- 
lated emission laser in which the three-level atoms are 
initially prepared partially in 50:50 probability to be in 
the upper and lower energy levels is considered (detailed 
description of the physical model can be found for in- 
stance in [Uj]). It is not difficult to envision at this junc- 
ture that the pumping radiation has associated phase 
fluctuations due to instability in the laser. It has been 
reported elsewhere that the effect of such fluctuations on 
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the generated entanglement is quit e perceivable specially 
when there is initial radiation |13l [171 [l9j . Despite this 
evidence, the phase fluctuations related to the imperfect 
preparation is considered here. 

In addition, in connection with recent observation that 
the rate at which the coherent superposition decays con- 
tributes significantly towards the evolution of the entan- 
glement [f9( and a proposal for the importance of includ- 
ing the rate of dephasing in the analysis , the rate at 
which the coherence superposition decays is assumed to 
be greater than the atomic decay rate. This assumption 
principally addresses the concern that the physical phe- 
nomena like vacuum fluctuations and broadening of the 
atomic levels can lead to loses of coherence without nec- 
essarily resulting atomic transitions. Taking these as mo- 
tivation, the effect of the processes that induce coherence 
(initial preparation and external pumping) in one hand 
and the processes that destroy coherence (dephasing and 
phase fluctuation) on the other hand on the evolution of 
the two-mode squeezing and intensity of the cavity radi- 
ation is investigated. 

In the study of the time evolution of similar system 
when many parameters are involved, the usual approach 
is numerically calculating the required correlations using 
either the characteristic function [j| [l3| or the rate equa- 
tion [TBI . Ht| . However, in this paper, the dynamics of 
the cavity radiation is analyzed following a straightfor- 
ward analytic approach. In the first place, the rate equa- 
tions corresponding to the expectation values of the radi- 
ation modes ((d) and (b)) are derived applying the master 
equation. Then the resulting expressions are converted 
into equations of c-number variables associated with the 
normal ordering. After that the obtained coupled dif- 
ferential equations are analytically solved employing the 
procedure outlined in [2(| and reference there in. Once 
the solutions are obtained various correlation are deter- 
mined with the aid of the correlations of the pertinent 
Langevin noise forces. 



II. EQUATIONS OF EVOLUTION 

Interaction of a pumped nondegenerate three-level cas- 
cade atom with a resonant two-mode cavity radiation can 
be described in the rotating-wave approximation and the 
interaction picture by the Hamiltonian of the form 

H = ig[a\a) (b\ - \b) (a\tf + b\b) (c| - |c) (b\tf] 

+ ij[\c)(a\-\a)(c\], (I) 

where O is a real-positive constant proportional to the 
amplitude of the driving radiation and g is a coupling 
constant chosen to be the same for both transitions, d 
and b are the annihilation operators that represent the 
two cavity modes. In the cascade configuration, the tran- 
sition from upper energy level \a) to the intermediate 
level 1 6) and from level \b) to the lower energy level \c) are 



taken to be resonant with the cavity radiation, whereas 
the transition \a) -H> \c) is electric dipole forbidden. Al- 
though the pumping laser undoubtedly has a bandwidth 
that consequently leads to phase fluctuations, this issue is 
not considered here to limit the involved rigor. However, 
due to emerging various quantum effects, it is assumed 
that the atoms can only be initially prepared in a partial 
maximum coherent superposition of the upper and lower 
energy levels. 

In view of these, the initial state of the three-level atom 
is taken to be 



|* A (0)> = —[\ a ) + e i *\c)], (2) 



where ip is an arbitrary phase randomly distributed about 
a fixed mean value ipa. As a result, the phase can be de- 
fined as ip = ipo + Sip where Sip is taken to be small 
random fluctuations around ipa which can be adjusted at 
will by proper choice of the phase of the cavity radiation 
;15|. Hence one can set tpo = and take ip as fluctuations 
about central value 0. In line with Eq. ([2]), the corre- 
sponding initial reduced atomic density operator appears 
to be 

p A (0) = ±[\a)(a\ + e-^\a)(c\ + e*\c)(a\ + \c)(c\] , (3) 

where e ±lip /2 stand for the initial atomic coherence. 

In setting up the laser, the atoms prepared in this man- 
ner are assumed to be injected into the cavity at a con- 
stant rate r a and removed after sometime which is long 
enough for the atoms to decay spontaneously to energy 
levels that are not involved in the process. The rate at 
which the atoms decay through spontaneous emission T 
is assumed to be less than the rate at which the atomic 
coherence superposition decays 7. Assuming a good cav- 
ity limit, the case in which the cavity damping rate k 
is much smaller than various channels of atomic decay 
rates (r and 7), the adiabatic approximation scheme is 
employed. For the sake of convenience, these damping 
constants are taken to be the same for different modes 
and energy levels as the case may be. 

Moreover, the linear analysis is applied in which all 
terms are limited to the square of the coupling constant 
g so that the resulting differential equations can be an- 
alytically solvable. In the study of the nondegenerate 
three-level laser, the linearization procedure is usually 
applied in connection to the fact that the quantum fea- 
tures are by large associated with the coherence induced 
in the cascading process. With this background, apply- 
ing the linear and adiabatic approximation schemes in 
the good cavity limit, the time evolution of the reduced 
density operator for the cavity radiation is found follow- 
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ing the procedure outlined in [12, 20] to be 
AC 



— = — [2apa ' — a ' ap — pa 1 a J 



25 



where 



\2a) pa — pad) — ad) p~\ 



i f ^ + «^ 26p&t _ pgtg _ b r bp 



AD + 

2B 

AD- 

2B 

AE+ 

2B 
AE_ 



2B 



2B 



bpw — a) pa — b'bp + da) p 
bptr — a) pa — pl)b + pad) 
a) plr — b'a)p + bpa — abp 
a) ply — ply a) + bpa — pab 

l)d)p — ply a) — abp + pab 



A = 



2r a g 2 



B=(4 + C 2 )(1 + C'C), 



C-2(C' 2 +x), 



(4) 



(5) 



(6) 



(7) 



mechanisms have competing effect in which the quantum 
features that is lost due to faulty preparation setups can 
be compensated by adjusting the amplitude of the exter- 
nal driving radiation. 

On the other hand, as clearly demonstrated elsewhere 
[H EH H3| i the stochastic differential equations associated 
with the normal ordering of the cavity mode variables are 
found making use of Eq. (0| to be 

d Bk + A[{2C + Qe- d -2(C + X)] m 
dt a{t) = 2B a{t) 



A [C'(l + CO - (2 - C'C)e- e ] 



2B 



(3*(t) + f a (t), 
(10) 



>^- ^ +A[(2c,+ r e+2(c,+x)] ^ 



2B 



A[c(i+(C) + (2-ec> 

2B 



-a*(t) + f b (t), 
(11) 



where f a (t) and fb(t) are the stochastic noise forces whose 
various correlations turn out to be 

(fa(t')f:(t)) = ^ [2C 2 + 2 X - (2C + Qe- 6 } 6{t - t'), 

(12) 



(Mt')fb(t))=0, 



(13) 



At = (2C + ()e ±llp , 



(2 - C'C)e 



(8) 



(9) 



in which £ = Cl/-y, (' — il/T, and \ = 7/r. The contri- 
bution of the cavity damping which corresponds to the 
coupling of the cavity modes with the two-mode vacuum 
reservoir via the coupler mirror is incorporated following 
the usual standard approach [2l], |22| ■ 

Comparison reveals that this master equation has a 
fundamental difference in form from the earlier report 
of similar scheme [llj]. It looks from the outset that 
the gain of mode a and the lose of mode b are equal. 
Even though the terms associated with C have the same 
form with earlier reports where such deduction is possi- 
ble, in Eq. some of these terms are unambiguously 
included in D± . As recently proposed, such modification 
in the master equation emanates from the changes in the 
subsequent atomic transitions due to phase fluctuations 
introduced during preparation [l2j ]. It is also vividly no- 
ticeable that the last term in the master equation closely 
related to the manifestation of the nonclassical features 
is solely attributed to changes of sign in E± due to phase 
fluctuations. Moreover, critical observation of the mas- 
ter equation reveals that the sign in front of b^aJp and 
pab in E± and the last term is different. This entails 
that the external driving and partial initial preparation 



(h(t')fa(t)) = ^ [C'(i + CO + (2 - C'Oe- 6 ] 5(t - t'), 

(14) 



(mt')fa(t)) - (fa(t')fa(t)) = (fb(t')f b (t)) = 0. (15) 

While deriving the above equations, a stochastic average 
is taken over the phase since it is assumed to fluctuate 
about the central mean value (set to for convenience). 
From practical point of view, addressing the contribu- 
tion of every phase change appears to be unrealistic. In 
connection to this, assuming that the phase undergoes 
Gaussian random process [lj, [TJ|, [Tj| , it would be more 
appropriate if the phase fluctuation is used instead of the 
actual phase. With this understanding and employing 
the fact that for Gaussian variables 



(exp±iS(p) = exp — (5tp 2 /2) 



(16) 



exp±i(/3 is replaced by exp— 9 [2J]. It is not difficult to 
realize that for Gaussian random process (6<j)) is zero, 
hence 6 = (8 2 (j>/2) represents the deviation which is gen- 
erally designated as phase fluctuation. 

In order to see the time evolution of each variable, 
these coupled differential equations are solved following 
the approach outlined in [12, [2(| . The solutions are 

a(t) = F + (t)a(0) + G+(t)p*(0) + H + (t) + J + (t), (17) 
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/3{t) = F-(t)/3{0) + G-(t)a*(0)+H-(t) + I-(t), (18) 
where 

F ± (t) = I [(1 ±p)e-»- t + (1 T P )e^ +t ] , (19) 



■ii-t] 



(20) 



H + {t) = \f\l+p)e 



M-(i-t') 



+ (1-P)« 



-M+(t-t')l fJ+>\J+' 



H.{t) = \j\l-p)t 



}fa(t')dt', (21) 



-/i-(i-t') 



+ (l+p)e-' J +( t -*')]/ 6 (t')^', (22) 
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/ !/) = y ^ [e-"+(*- t ') - e^*-*')]/* (*')< (23) 



J_(t) = q — J [ e -/^(*-*') - e -/*-(*-0]/*(f')dt' | (24) 
with 

M± = ^ + -4[(2C' + C)e- 9 



2 2B 



±\/C' 2 (l + CO 2 + 4(C' 2 + x) 2 - [(2 - C'C)e 



(25) 



2[C' 2 + x] 

VC' 2 (1 + CO 2 + 4(C' 2 + x) 2 - [(2 - C'C)e-T ' 



(26) 



<7± 



-C'(l + C'C)± [(2- C'C)e- e ] 



\/C' 2 (l + CO 2 + 4(C' 2 + x) 2 - [(2 - C'C)e~T 

(27) 



It perhaps worth mentioning that Eqs. (|17[) and (fT8|) 
along with the associated parameters are used to calcu- 
late various quantities of interest. 



III. TWO-MODE QUADRATURE SQUEEZING 

The two-mode squeezing for coherently pumped cor- 
related emission laser has been studied at steady state 
when the atoms are assumed to be perfectly prepared 
initially [llj]. Taking all atomic decay rates as equal, it 
has been reported that the two-mode squeezing exists 



either for smaller or larger amplitudes of the driving ra- 
diation. In connection to the restriction associated to 
steady state solutions, the analysis has been limited to 
smaller values of the linear gain coefficient. Comparison 
shows that the degree of squeezing is better for O = O.57. 
In this section, following a similar approach, the evolu- 
tion of the two-mode quadrature squeezing is evaluated 
and comparison with the steady state case is made when- 
ever possible. In addition, the effects of phase fluctuation 
and dephasing on the degree of squeezing in relation to 
the pumping mechanism is investigated. The cases for 
which these phenomena, believed to degrade nonclassical 
features, aid in enhancing the squeezing is also looked 
for. 

Generally, a two-mode cavity radiation can be de- 
scribed by an operator 



1 



(28) 



where a and b represent the separate modes. With this 
consideration, the squeezing properties of the cavity ra- 
diation can be studied applying the quadrature operators 
defined by 



c + = C T + c 



and 



(29) 



(30) 



where the corresponding variance can be readily obtained 
in terms of c-number variables associated with the normal 
ordering as 



Ac| = l + i[2(a*a)±(a* 2 )± 



a 2 ) - (a*) 2 - (a) 2 



± (/T 2 ) ± ((3 2 ) + 2(/T/3) - (/T) 2 - (f3) 2 ] ± (a*/?*) 

+ (a*(3) + (a/3*) ± (a/3) T (a*) (a) - («*><£*) 

T (a*)((3) T (a)(n - (a)((3) T (W)- (31) 



The various correlations in Eq. pip can be determined 
using Eqs. (TT7|) and (|T5)) . In this line, assuming the 
cavity to be initially in a two-mode vacuum state and 
the noise force at time t does not statistically related to 
the cavity mode variables at earlier times, one can readily 
verify that 



(a*(t)a(t))=A 
+ A 



(«(*)> = (P(t)) = 0, 
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FIG. 1: Plots of the quadrature variance (Ac?.) of the cavity 
radiation for n = 0.5, 7 = V, 9 — 0, Q, = O.57, and different 
values of A. 



(a 2 (t)) = ((3 2 (t)) = (a*(t)f3(t))=0: (36) 

where 

L = 2(' 2 +2x-(2C + ()e- 9 , (37) 

M = C'(l + C'C) + (2-C'C)e- e . (38) 

Hence on account of Eqs. (1321) and (|36|) . the quadrature 
variances of the cavity radiation reduces to 

Ac 2 ± = l + (a*(t)a(t)) + (P*(t)l3(t)) 

±[(a*(t)l3*(t)) + (a(t)(3(t))}- (39) 

In order to understand the dependence of the evolution 
of the two-mode squeezing on the linear gain coefficient, 
amplitude of the external driving radiation, the rate 
of dephasing, and the phase fluctuation in depth, case 
by case approach is preferred. Throughout this work, 
two-dimensional plots are chosen instead of the three- 
dimensional surfaces for the sake of clarity. It is worth 
noting that the main aim is geared to show the nature 
of the evolution so that comparison can be easily made 
notwithstanding the temptation that the maximum ob- 
tainable degree of squeezing can be optimized by adjust- 
ing the involved parameters. 

A. Dynamics of the two-mode squeezing 

Let us first take the case when there is no phase fluc- 
tuation {9 — 0) and the atomic damping rate via every 
channel is the same (7 = F). In the following, the plots 
of the variance of the quadrature operators versus time 
is presented using Eqs. (JM]),®,®,©,®,©, 

(EB, 42$, and d3U). 
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FIG. 2: Plots of the quadrature variance (Ac?.) of the cavity 
radiation for n = 0.5, 7 = V, 9 — 0, Q = 2.57, and different 
values of A. 

It is clearly shown in Figs. HJ HI and |3] that the cavity 
radiation exhibits two-mode squeezing starting from the 
inception of the injection of the atoms. It is not difficult 
to notice that, in the early stages of the operation, the de- 
gree of squeezing rapidly increases with time irrespective 
of the values of the linear gain coefficient and amplitude 
of the external pumping radiation. This can be viewed 
as the injection of more atoms into the cavity leads to 
the generation of more correlated photons. This result 
in other words indicates that in the early stages of the 
amplification, the nonclassical properties of the radiation 
rapidly increase with time. However, as time passes by 
the thermal fluctuations arising due to heating begins to 
consume the established nonclassical features and hence 
the squeezing starts to wear away. 

The values of fi/7 are chosen based on the earlier stud- 
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FIG. 3: Plots of the quadrature variance (Ac?_) of the cavity 
radiation for k, = 0.5, 7 = V, 9 = 0, Q, — IO7, and different 
values of A. 



ies where the two-mode squeezing for similar case is found 
to be strong at d — O.57, very poor when S7 — 2.57, an d 
some what intermediate when Q = IO7 at steady state 
pjjj]. In light of this, critical comparison of results shown 
in Figs. [U [2j and [3] asserts that the degree of squeezing 
is better for fi = O.57 in the early stages of the evolu- 
tion. Quiet remarkably, contrary to steady state case, it 
is possible to witness squeezing for £7 = 2.57, although 
the degree of squeezing is still small when compared to 
the other two. Extrapolating this representative result 
evinces that it is possible to obtain squeezed light from 
this system for all values of the amplitudes of the external 
driving radiation in the early stages of the amplification 
process. It is also worth noting that the time at which 
a maximum squeezing can be manifested depends on the 
linear gain coefficient and the amplitude of the driving 
radiation. Critical scrutiny of Figs. [TJ [5J and [3] reveals 
that for similar values of A, the time at which the degree 
of squeezing would be maximum is the smallest when 
fl = 2.57 and largest when SI ~ IO7. Most probably this 
can be related to the absence of squeezing when SI = 2.57 
at steady state. The other issue that may worth mention- 
ing is the enhancement of the degree of two-mode squeez- 
ing with the linear gain coefficient at early stage of the 
operation and then the situation changes after maximum 
squeezing has been reached. This result also suports the 
explanation provided in relation to the thermal fluctua- 
tions earlier. 



B. Effect of the phase fluctuation on the dynamics 
of the two-mode squeezing 

The dependence of the evolution of the two-mode 
squeezing on the phase fluctuation for different values 
of the amplitude of the external driving radiation and 
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FIG. 4: Plots of the quadrature variance (Ac 2 .) of the cavity 
radiation for k — 0.5, 7 = T, A — 10, = O.57, and different 
values of 9. 
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FIG. 5: Plots of the quadrature variance (Ac?,) of the cavity 
radiation for k — 0.5, 7 = T, A — 10, = 2.57, an d different 
values of 9. 

It is not difficult to see from Figs. 0] [5j and |H] that 
the two-mode squeezing increases with time in the early 
stages of the operation and wears away as time progresses 
for the reasons already explained. It is clearly shown in 
Figs. S]and[5]that the degree of squeezing decreases with 
the enhancement of the phase fluctuation. Contrary to 
this, in Fig. [51 the degree of squeezing is found to increase 
with the phase fluctuation. The same phenomenon has 
been repeatedly observed when ft is taken to be much 
greater than 7. Critical comparison of the results shown 
in Figs. |3] and [5J for A = 10, indicates that the two- 
mode squeezing is better when the atoms are initially 
prepared in a partial coherent superposition. This result 
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FIG. 6: Plots of the quadrature variance (Ac?.) of the cavity 
radiation for k = 0.5, 7 = T, A = 10, fi = IO7, and different 
values of 9. 



FIG. 7: Plots of the quadrature variance (Ac?.) of the cavity 
radiation for n — 0.5, 9 = 0, A = 10, fi = O.57, and different 
values of 7/T. 



entails that the squeezing that could have been lost due 
to phase fluctuation can be compensated for better by 
the external driving mechanism and vice versa via creat- 
ing indirect roots for spontaneous atomic emission. The 
other noteworthy aspect of the phase fluctuation, in the 
strong driving limit, is the relaxation of the time at which 
the maximum squeezing is manifested and the time for 
which the radiation remains in the squeezed state which 
might be useful in practical situation in connection with 
the smallness of the actual time. 



C. Effect of dephasing on the dynamics of the 
two-mode squeezing 

In the following, upon fixing A — 10 and varying the 
values of 7/T the effect of dephasing on the evolution 
of the two-mode squeezing is analyzed. Here the rate 
at which the coherent superposition decays is assumed 
to be greater than the other atomic decay rates. This 
assumption appears to be more realistic than assuming 
the two variant of atomic dampings as equal, since the 
induced atomic coherent superposition can be lost via 
quantum processes like vacuum fluctuations and atomic 
broadening without directly involving atomic transitions. 
Basically all phenomena related to atomic broadening, 
whereby the coherence between the lower and upper en- 
ergy levels can be lost without necessarily one of them 
decaying to other levels, can lead to enhanced dephasing. 
Apparently, the effects of dephasing on the evolution of 
the two-mode squeezing is believed to be an interesting 
issue by its own right. 

It is possible to infer from Fig. [7] that the degree of 
two-mode squeezing decreases with increasing rate of de- 
phasing for smaller amplitudes of the driving radiation. 
Similar to earlier cases the degree of two-mode squeezing 
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FIG. 8: Plots of the quadrature variance (Ac?.) of the cavity 
radiation for n — 0.5, 9 = 0, A = 10, = IO7, and different 
values of 7/r. 



increases with time at earlier stages and then decreases 
with time. It is worth to note that in addition to the 
squeezing the time at which the maximum squeezing ob- 
tained is also decreases with dephasing. Therefore, the 
practical advantage of dephasing in weak driving limit 
is far from encouraging. In the case that the phase can 
be locked via for instance coherent population trapping 
through two-step external pumping, that means when 
7 < r, the two- mode squeezing significantly increases 
with decreasing 7/r. A similar result has been recently 
reported for quadrature entanglement [Iff. Contrary to 
this, as clearly shown in Figs. [8] and HI the degree of 
two-mode squeezing increases with the rate of dephasing 
when the atoms are externally pumped with a strong ra- 
diation. It is straightforward to see that in Fig. [5] the 
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FIG. 9: Plots of the quadrature variance (Ac?.) of the cavity 



radiation for k 
values of 7/r. 



0.5, t = 0.85, A = 10, fi = IO7, and different 



phase fluctuation is varied when 7 = T, but in Fig. [8]thc 
dephasing is varied when 9 = 0. In these cases, the degree 
of squeezing is found to increase with phase fluctuation 
and dephasing. On the other hand, as clearly shown in 
Fig. |9l both dephasing and phase fluctuation are made to 
vary. It turns out that the degree of two-mode squeezing 
increases with both. This result testifies that the assump- 
tions that have been made for limiting the mathematical 
rigor have overshadowed the obtainable nonclassicality of 
the radiation. What surprises most is these assumptions 
have been thought to increase the same which of course 
is true in weak driving limit or when there is no driving 
at all. 



IV. MEAN NUMBER OF PHOTON PAIRS 

In relation to the brightness of the squeezed light to 
be generated, it is necessary to evaluate the mean pho- 
ton number which is directly associated with the inten- 
sity of the radiation. In order to attest the potential of 
the system under consideration as a reliable source of a 
two-mode squeezing, the evolution of the photon number 
for the cavity radiation would be studied. Genrally, the 
mean number of photon pairs describing the two-mode 
cavity radiation can be defined as 



N 



f(t)c(t)), 



(40) 



where c(i) is the annihilation operator given by (|28l) . 

On the basis that the operators in Eq. (|4U|) are already 
put in the normal order, it is possible to rewrite it in 
terms of c-number variables associated with the normal 
ordering as 



1 



Since (a*(t)a(t)) and (f3* (t)f3(t)) represent the mean 
photon number in mode a and mode b, N is interpreted 
as mean number of photon pairs. The actual intensity 
of the radiation in the cavity can be taken as twice of N 
without essential lose of generality. In the following, to 
save time and space, only certain representative points 
are selected. The parameters in the analysis are fixed 
to be the same as the corresponding two-mode squeezing 
except in certain cases where some changes are made for 
the purpose of clarity. 
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FIG. 10: Plots of the mean number of photon pairs (N) of 
the cavity radiation for k = 0.5, 7 = F, 8 — 0, fl = O.57, and 
different values of A. 




N=-[(a^t)a(t)) + (P*(t)0(t))]- 



(41) 



FIG. 11: Plots of the mean number of photon pairs (N) of 
the cavity radiation for k = 0.5, 7 = V, A — 10, = IO7, and 
different values of 8. 

It is not difficult to see from Figs. [TU1[TT1 and mi that 
the intensity of the radiation increases with time irrespec- 
tive of the changes in the parameters under consideration. 
Unlike the two-mode squeezing, this trend prevails even 
at larger values of time. This outcome mainly evinces 
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FIG. 12: Plots of the mean number of photon pairs (N) of 
the cavity radiation for k — 0.5, t = 0.85, A — 10, = IO7, 
and different values of 7/r. 



that the more the atoms are injected into the cavity, the 
more photons are available in the cavity. The thermal 
fluctuations arising from heating of the cavity that has 
been found to consume the nonclassicality has no effect 
on the number of the photons to be generated. In line 
with this, it is possible to envision that this system can 
generate sufficiently bright light at steady state as al- 
ready predicted in earlier communications fl8| . It is also 
vividly presented in Fig. [10] that the intensity of the radi- 
ation increases with linear gain coefficient in low driving 
limit (Q < 7). Since the atoms are initially prepared 
partially to be 50% in the upper energy level, the more 
atoms are injected into the cavity implies more radiation 
would be emitted when the amplitude of the external ra- 
diation is too small to reverse the trend. Moreover, as 
clearly shown in Figs. [TT1and[T2l the intensity of the ra- 
diation increases with the rate of dephasing and degree 
of phase fluctuation for selected parameters. 



V. DISCUSSION AND CONCLUSION 

In this contribution, analysis of the time evolution of 
the two-mode squeezing and intensity of the cavity ra- 
diation of a coherently pumped two-photon correlated 
emission laser is presented. The effect of dephasing re- 
sulting from the involved quantum phenomena like vac- 
uum fluctuations and atomic broadening along with the 
phase fluctuation corresponding to incapability of prepar- 
ing the atoms initially in a perfect 50:50 probability to 
be in the upper and lower energy levels on the dynamics 
of the squeezing are the main issues. In order to see the 
contribution of the driving radiation in depth, three rep- 
resentative cases are selected: the regime of weak driving 
limit (ft = O.57), the case where the squeezing is expected 
to be minimum (O = 2.57), and a strong driving limit 
(fi = IO7). In order to enhance clarity of exposition, the 



presentation is divided into three main categories. In the 
first place, the evolution of the two-mode squeezing due 
to changes in the linear gain coefficient and amplitude of 
the driving radiation is treated. Secondly, effect of the 
phase fluctuation along with the variations in the ampli- 
tude of the external driving radiation on the evolution 
of the two-mode squeezing and intensity of the radiation 
is investigated. Lastly, effect of dephasing in relation to 
modifications in amplitude and phase fluctuation on the 
same is addressed. Although this study is limited to only 
few choices of the involved parameters, the obtained out- 
comes entail that the approach adopted can lead to a 
more tractable results than the approximated numerical 
analysis can offer. 

It turns out that the degree of two-mode squeezing 
and the intensity of the generated radiation increase with 
time at the early stages of the amplification irrespective 
of the alterations in the parameters under consideration. 
However, this tendency changes as time progresses spe- 
cially in case of the squeezing. It may worth mentioning 
that a similar conclusion has been drawn earlier under 
various conditions [f| [l3l - [l5| . On the other hand, the 
degree of two-mode squeezing begins to rapidly decrease 
with time whereas the intensity of the radiation remains 
to increase in later stages of the operation. It is not 
difficult to envision that with increment of time more 
atoms under went spontaneous emission while traversing 
across the cavity. Undoubtedly, the more there are emit- 
ted photons available in the cavity, the more would be 
the temprature of the cavity. Due to this heating, the 
atoms on the wall of the cavity starts to vigorously vi- 
brate which causes additional thermal fluctuations. It 
is this thermal fluctuation that consumes the correlation 
to be established between the modes and consequently 
leads to inhibition of the nonclassical features. Since the 
thermal fluctuations do not have a direct impact on the 
number of emitted photons, the intensity of the radiation 
continues increasing with time as long as the atoms are 
injected into the cavity. Even though it is not evident in 
this study, it has been reported elsewhere that the gener- 
ated radiation regains its nonclassical behavior at steady 
state since the coherence due to correlation of emitted 
photons slowly overcomes the effect of the thermal fluc- 
tuation through time [l|, [l8[ . 

The dependence of the degree of two-mode squeezing 
on the dephasing and phase fluctuation is found to be dif- 
ferent in the weak and strong driving limits. In the weak 
driving limit, the two- mode squeezing decreases with 
both as usually expected. However, in the strong driving 
limit, the two-mode squeezing increases with both de- 
phasing and phase fluctuation. Evaluating the situations 
from the root betoken that, in the strong driving limit, 
the coherence induced by external radiation is too strong 
to allow spontaneous emission in the cascading process. 
The evidence may be traced to a significant decrement 
in the mean photon number with the amplitude of the 
driving radiation in a similar case at steady state @. 
Hence dephasing and phase fluctuation characteristically 
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weaken this coherence whereby provide an indirect elec- 
tron path way which enhances the emission of correlated 
photons. 

Based on the assumption that there is no interaction 
between the injected atoms, that is clearly manifested in 
Eq. (f36|) . the relation between the two-mode squeezing 
and the entanglement measure due to the criterion fol- 
lowing from Duan et al. [25| for general nondegenerate 
three-level laser has been established [f| [23[ . In view of 
these reports, the results obtained for the degree of two- 
mode squeezing are equally applicable for the quadrature 
entanglement. With this understanding, it is possible 
to conclude that the technical inability of preparing the 
atoms in a perfect coherent superposition and the fail- 
ure to control the destruction of coherence by quantum 
fluctuations provide a means for enhancing the quantum 
features including entanglement. If such unexpected out- 
comes can be properly harnessed, hopefully they support 
the pursuit for generating strong light with robust non- 
classical features. Moreover, the possibility of increasing 
the time at which maximum squeezing is observed and 
the time in which the light remains in the squeezed state 



by manipulating the phase fluctuation can be utilized 
in the practical realization of this system. Furthermore, 
the increment of the intensity of the radiation with time 
complements the advantages of this system as a poten- 
tial candidate for the generation of the light with strong 
nonclassical features. Since the assumptions that have 
been made for some particular interest (which could have 
posed a sever restriction) are lifted, it is hoped that the 
system under consideration can be tested in a more re- 
laxed manner with anticipated far better outcomes. With 
no doubt, this would be a positive signature in the pursuit 
of generating strong nonclassical light from this system. 
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